UICHEP-TH/96-18 
quant-ph/96TTmU 



Accuracy of Semiclassical Methods 
for Shape Invariant Potentials 



Marina Hruskaft, Wai-Yee KeungQ and Uday SukhatmeQ 

Department of Physics, University of Illinois at Chicago, 
845 W. Taylor Street, Chicago, Illinois 60607-7059. 



Abstract 

We study the accuracy of several alternative semiclassical methods by computing analyt- 
ically the energy levels for many large classes of exactly solvable shape invariant potentials. 
For these potentials, the ground state energies computed via the WKB method typically devi- 
ate from the exact results by about 10%, a recently suggested modification using nonintegral 
Maslov indices is substantially better, and the supersymmetric WKB quantization method 
gives exact answers for all energy levels. 
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1. Introduction 



A variety of semiclassical methods have been proposed and used for determining the energy 
levels of one-dimensional potentials. The standard WKB method is discussed in most quantum 
mechanical texts [|T], [2j but there are several recently suggested modifications based on supersym- 
metric quantum mechanics 0, energy dependent phase losses at the classical turning points with 
nonintegral Maslov indices f|, and other related approaches 0. As expected, all semiclassical 
methods yield good energy eigenvalues E n for large values of the quantum number n. However, 
their accuracy varies quite substantially for small values of n depending on the choice of potential. 
Also, many potentials only have a small number of bound states, so that the possibility of consid- 
ering large values of n does not exist. In this paper, we make a stringent test of the accuracy of 
various semiclassical methods by computing the eigenenergies of the ground state and other low 
lying states for many large classes of potentials which have the property of shape invariance || 
under supersymmetry transformations. We have chosen shape invariant potentials since (i) they 
are exactly solvable and all eigenvalues are explicitly known; (ii) the integrals appearing in the 
semiclassical quantization conditions can all be performed analytically; (iii) the nonintegral Maslov 
indices used in the recent semiclassical approach M of Friedrich and Trost (FT) can be expressed 
in terms of superpotentials. Our plan is to first review three semiclassical approaches and the main 
ideas involving shape invariant potentials. We will then compute the quantization condition inte- 
grals analytically and determine energy eigenvalues. The results are presented in Tables 1 and 2, 
and we give some concluding remarks on the accuracy of various semiclassical approaches. 

2. Semiclassical Quantization Conditions 

(i) WKB Quantization: The usual form of the semiclassical energy quantization condition 
(in units of % = 2m = 1) is |[ g 

f XR dxy/E - V{x) = (n + , /i = , (n = 0, 1, 2, ...) , (1) 

Jx L v 4 TT/2 

where the classical turning points xl and xr are given by V(xl) = V(xr) = E. The Maslov index 
/j, denotes the total phase loss during one period in units of tt/2. It contains contributions from the 
phase losses 4>l and 4>r due to reflections at the left and right classical turning points xl and xr 
respectively. In the standard WKB approach, one takes <pL = <pR = tt/2, and an integer Maslov 
index fi = 2 for all energy levels. This gives the familiar result (n + \)n for the right hand side of 
the WKB quantization condition. 

(ii) SWKB Quantization: Another semiclassical approach || which has been widely studied 
in recent years is based on the ideas of supersymmetric quantum mechanics 0. Here, the super- 
symmetric partner potentials V-(x) and V + (x) are given in terms of the superpotential W(x) by 
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V± = W 2 (x) ± W'{x). For the case of unbroken supersymmetry, VL(x) and V + (x) have degenerate 
energy levels except that has an additional level at Eq = 0. The corresponding ground 

state wave function ipQ \x) is related to the superpotential W(x) via 

W(x) = -^^ ; 4~\x)<x e-r^') . (2) 
% >(x) 

The supersymmetric WKB (SWKB) approach |j| results from combining the ideas of supersymme- 
try with the lowest order WKB method. The SWKB quantization condition is 

f X " dx^jE(-) - W 2 {x) = mr , (n = 0, 1, 2, ...) , (3) 

Jx' L 

where the two turning points x' L and x' R are given by W(x) = ±y/E(~). Note that for n = 0, the 
turning points coincide and the SWKB quantization condition gives the exact result Eq~^ = for 
the ground state energy. 

(iii) Friedrich-Trost Quantization: This very recent proposal Q makes use of the stan- 
dard quantization condition [eq. (1)] with nonintegral, energy-dependent Maslov indices \i. More 
specifically, the phase loss is taken to be given by 

(f) L ip'(x L ) (j) R ip'(x R ) 

tanl — ) = ——. — r , tan( — ) = — - , (4) 

where k = \JE — V min and V min is the minimum value of the potential V(x). In Ref. Q], it was 
suggested that one could use the lowest order WKB wave function in eq. @ in order to determine 
the phase losses 4>l,4>r for any practical application. Indeed, it was shown that for power law 
potentials x p with p = 4, 5, 6, this method gave better numerical results for the ground state 
energies than the standard WKB method, and also substantially improved wave functions. 

The FT approach of using the WKB wave function in eq. (Q) is rather cumbersome. It can be 
significantly simplified for the ground state by using eq. (|2|) • The phase losses (pL, 4>R can then be 
re-written as 

, Al, W{X L ) (f)R W(x R ) 

tan(y) = ^— , tan(— ) = — - — . (5) 

In this paper, we will assume eq. (5) to be valid for all energy levels in computing eigenenergies by 
the FT approach. 

3. Shape Invariant Potentials 

Given a superpotential W(x, a Q ) depending on a set of parameters a , the supersymmetric 
partner potentials V±(x, a ) are given by 

V±{x,a ) = W 2 {x,a ) ±W'{x,a ) . (6) 
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These partner potentials are shape invariant if they both have the same x-dependence Upton a 
change of parameters a\ = /(ao) and an additive constant R(ao). The shape invariance condition 
is 

V+(x,a ) = V-(x, ai) + R(a ) . (7) 

This special property permits an immediate analytic determination of energy eigenvalues and 
eigenfunctions 0. For unbroken supersymmetry, the eigenvalues are 

n-l 

E { - ] = 0, £H = £ R ( a k) ■ (8) 

fc=0 

Many aspects of the bound states and scattering matrices of shape invariant potentials have been 



studied including several choices |I| for the change of parameters ai = /(ao)- In this paper, 
we confine our attention to shape invariant potentials corresponding to a translational change of 
parameters. 

4. Computation of Energy Eigenvalues 

All known families of shape invariant potentials in which the change of parameters is a translation 
a\ = a + j3 are listed in Table 1. Names of the potentials and the corresponding superpotentials are 
given. Also tabulated is the minimum value V- m %n of the potential V-(x) and the position x min of 
the minimum. Subtracting V- m i n from V-{x) yields the tabulated potential V(x), whose minimum 
value is clearly V m i n = 0. The exact energy eigenvalues E% cact for V(x) coming from eq. (^) are 
also given. 

To assess the accuracy of various semiclassical approaches, the first step is to evaluate the two 
types of integrals appearing in the quantization conditions. We denote the integral in the WKB 
condition eq. (1) by I WKB and the integral in the SWKB condition eq. (3) by J SWKB _ The integrals 
can be handled analytically using contour integration in the complex plane taking special care of 
the singularities at infinity and the cut going between the turning points. The results are given in 
Table 1 in terms of the following expressions: 

h(a,b)=[ dyJ (y -a)(b-y) = J(6 - a) 2 ; 

J a O 

I 2 (a, b) = — J (y - a)(b - y) = -(a + b) - nVab , (0 < a < b) ; 
Jay z 



la y 

(a, b) = f *L-y/(y- a )(b-y) = -^=lVTT^VTT¥ -ab+ I] 1 ' 2 - n 
J a y + l v 2 



h(a,b) = I" ^— 2 ^{y - a) (6 -y) = ^[2-^(1 - a)(l - b)- sj(l + a)(l + b)] , (-K a < b < 1) ; 
Ja i — y z 
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h(a,b) = f - - y) = + !)(&+!)- " W " l ) " 2 ] >(!<«<*) 5 

J« y — 1 Z 

In all the above integrals, the limits a, b are real numbers with a < b. We have given explicit 
expressions for the above integrals since they are not easily available in standard integration tables. 
Once I WKB and J SWKB have been computed, one can apply the quantization conditions to see 
how accurate the WKB, SWKB and FT approaches are. For the WKB and SWKB approaches, 
it is possible to get complete analytic results - these are shown in Table 1. The FT approach is 
also mostly analytical, but the final computations need numerical work. Results corresponding to 
specific numerical choices of the parameters appearing in the potentials are shown in Table 2. 

As an illustrative example, consider the Rosen-Morse II (hyperbolic) potential, for which the 
superpotential is 

W(x) = At&nhax + — . 
With a change of variables y = tanh ax, the SWKB integral is 



1 - aJy' L l-y 2 V V A 2+{ A 2 V J 



with turning points given by Ay' + ;f = ±V E(~\ One then sees that the integral is -^(yi,?/^)- 
Substitution into the SWKB quantization condition l SWKB = nn and solving for E^ gives 

E i-)SWKB { 2 { A . na) 2 . I_JL & 



A 2 (A - no) 1 

which is the exact answer for all energy levels ! Similar steps give the WKB integral to be I WKB = 
y/ 'A(A+a) j^y L ^ y R ^ w h er g the turning point are given by 

Substitution into the WKB quantization condition I WKB — (n + i)7r and solving for the energy 
gives 



E™ KB = A(A + a) - (^A(A + a) - ^ - no) + 



B 2 B 2 



A ( A + a ) (yjA(A + a) - % - no) 



The WKB approach does not give the exact eigenvalues. The full energy computation for the FT 
quantization condition is harder to carry out analytically. For the numerical choice a = 1, A = 
2, B — 1, we see from Table 2 that the ground state energy Eq T is lower than the exact energy by 
3.2% whereas E^ KB is higher than the exact energy by 9.7%. 

5. Conclusion 
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We have given a complete analytic treatment of the energy levels of shape invariant potentials 
(with a translational change of parameters) for various semiclassical quantization conditions. As 
expected from previous work 0, the SWKB energy levels are exact. The WKB energy levels 
are exact for the one-dimensional harmonic oscillator and the Morse potentials only. For other 
potentials the results are not exact, and this has historically led to ad hoc Langer corrections JTT 



Typically, one sees from Table 1 that the ground state energy from the WKB method deviates from 
the exact result by about 10%. The new semiclassical approach of Friedrich and Trost is exact only 
for the one-dimensional harmonic oscillator and no other shape invariant potential. However, in 
general, it gives significantly more accurate ground state energies than the WKB method. 

This work was supported in part by the U. S. Department of Energy. 



Table Captions 

Table 1: List of all shape invariant potentials and their eigenvalues. Analytic expressions for 
the integrals in the WKB and SWKB quantization condition are given, along with the energy 
eigenvalues. For these potentials, the SWKB results are always exact, whereas the WKB results 
are exact only for the harmonic oscillator and Morse potentials. 

Table 2: Comparison of the exact ground state energies of shape invariant potentials with results 
from the WKB and Friedrich- Trost method. The percent errors are also shown. The SWKB results 
are not shown, since they are always exact for the potentials under consideration. 
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Name 
of 

potential 


Super 
potential 

W(x) 


V_„„„, the 
minimum of 

v_ = w 2 - w 


— V—min 


v(x) = v_ (i) 

-V_ mi „, Note: 
V(x min ) = 


gexact j^SWKB 


jWKB _ 

QdxjE-V(x) 


jSWKB _ 

QdxjE-W 2 (x) 


j^WKB 


Shifted oscillator 


-LOX 


— -111 





V* 2 


(n+\)u> 


irE/io 


TlE/lH 


(n + \)u> 


Three 
dimensional 
oscillator 


1 (1+1) 


-(' + !)* 


[4l(l + l)/ui 2 ]i 


\uj 2 r 2 + l{l + l)/r 2 

-10^/1(1 + 1) 


2mj - V- m i n 




7TE 

2w 


2nw + w 


Coulomb 


e' _ (1+1) 
2(^+1) r 


e" 

il(l-\-\) 2 


2l(l + l)/e 2 


e' , 1(1+1) 

r r 2 

4 

+ 4(P + I) 


<L< 1 

(„+i+i) j ) 


■Ke 2 /^/e 4 /(l 2 + l)-4E 


■Ke 2 /^/e 4 /(l + l) 2 -4E 
-tt(Z + 1) 


iKP + i)- 1 

-(r^ + ^/P ^ +7+i)- 2 ] 


Morse 


A 

-Be- ax 


->la - a 2 /4 


I l0 S A+a/2 


(A + a/2 - Be- ax ) 2 


(.4 + a/2) 2 
-(^4 - na) 2 




1(A-JA 2 -E) 


(A + qy 
-(A - na) 2 


Rosen- 
Morse II 
(hyperbolic) 


A tanh ax 

+B/A 
(B < A 2 ) 


B 2 /A 2 - Aa 
-B 2 /(A 2 + Aa) 


tanh" 1 


B 2 /(A 2 + Aa) 
+A(A + a) tanh ax 
+-2_Btanho;^ 


A 2 + Aa - (A - na) 2 

. B 2 B 2 
~ I ~A 2 + Aa (A-na) 2 


y/A(A + a)I 4 (y L ,y a )/a 
y L + y R = -2B/(A 2 + Aa) 

B 2 E 
VLVR — (A 2 + An) 2 A 2 +An 


Ah(y' L ,y' R )/a 
y' L + y' R = -2B/A 2 
v' v' — SL S. 

ULUR — A' A 2 


a(l + 2n)y/A(A + a) 
-a 2 (i\ + n) 2 + B 2 /(A 2 + Aa) 
-B 2 /{VA 2 + Aa-%-na) 2 


Eckart 


— ^4 coth or 

+B/A 
( R ^ 4 2 "l 

(-D > A ) 


B 2 /A 2 + Aa 
-B 2 /(A 2 -Aa) 


1 coth" 1 ,, B , 

ct A z —Aa 


B 2 /(A 2 -Aa) 
+A(A - a) coth 2 ar 
— 2B coth ar 


A 2 - Aa - (A + na) 2 

, B 2 B 2 
~ I ~A 2 -Aa (A+na) 2 


^A(A- a)h{y L ,y a )/a 
Vl + Vr = 2B/(A 2 - Aa) 

— B 2 E 
VLVR — { A 2 -Aa) 2 A 2 -Aa 


Ah(y' L ,y' R )/a 

y' L + y' R = 2B/A 2 

,,l „l - B 2 E 
UlUr — A* A 2 


-a{l + 2n)y/A(A - a) 
-a 2 (\ + n) 2 + B 2 l(A 2 -Aa) 

R 2 /(</ A2 ^^-L^-L n /~v \ 2 
— lj J y\j i\ — s\OL -\- ~r TKXj 


Rosen- 
Morse I 
(trigonometric) 


— Acoi, ax 
-B/A 
(0 < ax < 7r) 


B 2 /A 2 - Aa 
-B 2 /(A 2 - Aa) 


- i cot" 1 ,, B , 


B 2 /(A 2 - Aa) 
+A(A - a) cot 2 ax 
+2B cot ax 


-A 2 + Aa + (A + na) 2 

+ A 2 -Aa (A+na) 2 


^/A{A-a)I 3 {y L ,y R )/a 
Vl + Vr = -2B/(A 2 - Aa) 

VLVR — (A 2 -Aa) 2 A 2 -Aa 


Ah(y' L ,y' R )/a 

y'L + y'R = -2B/A 2 

i/LilR — A 1 A 2 


+a(\ + 2n)^/A(A-a) 
+a 2 (| + n) 2 + B 2 /(A 2 - Aa) 
-B 2 /{VA 2 - Aa + f + na) 2 


Scarf II 
(hyperbolic) 


A tanh ax 
+-Bsech ax 
D = 2AB + aB 


A 2 -L 
C = B 2 — A 2 — Aa 
1L = -JC 2 + D 2 - C 




L + Csech 2 aa; 
+£>sech ax tanh ax 


L-(A- na) 2 


1{4Z--JL-E) 


±(A-^A 2 -E) 


L - (\/L - na - f ) 2 


Generalized 
Poschl- 
Teller 


A coth ar 
— Bcsch ar 
(A < B) 


A 2 -L 
C = A 2 + B 2 + Aa 
2L = C- VC 2 - D 2 


i cosh -1 ^ 
D = 2AB + «B 


L + Ccsch 2 or 
— Dcsch ar coth ar 


L-(A- na) 2 


l(VZ-s/L-E) 


l(A-^A 2 -E) 


L-(VL-na-f) 2 


Scarf I 
(trigonometric) 


A tan ax 
—B sec ax 
(A > B) 


-A 2 + L 
C = A 2 + B 2 -Aa 
2L = C+ -JC 2 - D 2 


1 s ; n -i _d 

a blU 2L 

D = 2AB - aB 


—L + Csec 2 aa; 
— -Dsec ax tan ax 

(M < f ) 


(A + na) 2 - L 


K-jl + e-Vl) 


1(JA 2 + E-A) 


(VZ + na + ^) 2 -L 



Table 1: List of all shape invariant potentials and their eigenvalues. Analytic 
expressions for the integrals in the WKB and SWKB quantization condition are 
given, along with the energy eigenvalues. For these potentials, the SWKB results 
are always exact, whereas the WKB results are exact only for the harmonic 
oscillator and Morse potentials. 



Name of 
Potential 


Choice of 
Parameters 


Super 
Potential 


Potential 

V(x) 


j^exact 
^j^exact^ 


jpWKB 
^0 


Percent 
error 


Fried rich. 
-Trost 

et 


Percent 
error 


oscillator 




±wa; 


jL0 2 X 2 


(1.5w) 




n 
u 




n 


± nree 
dimensional 
oscillator 


LO — L 
I = 1 


r 


(v^/r-r) 2 


9 1 71 
Z.l / ID 

(6.1716) 


o 

z 


o.U ( /o 




U.Dl /o 


Coulomb 


c — 1 

I = 1 


I _ 2 

4 r 


2(i-|) 2 


U.UDZO 

(0.09722) 


U.OD / / 


Q 1 7% 
y . 1 ( /o 


U.UOOoU 


i n 7% 

1U. / A) 


Morse 


py O 

A = 1 
B = 1 


1 - e- 2 * 


(2_ e -2a)2 


3 


3 





2.7587 


8.04% 


Rosen- 
Morse II 


py 1 

A = 2 
B = 1 


2 tanh x + \ 


i(l + 6 tanh a;) 2 


1.9167 


2.1030 


9.72% 


1.8549 


3.22% 


Eckart 


OL — 1 

A = 2 
B = 6 


3 - 2 cothr 


2 (cothr - 3) 2 


7 


6.511 


7.00% 


6.1540 


12.1% 


Tit riQOTl - 

Morse I 
(trigonometric) 


a = 1 
A = 2 
B = 2 


— 1 — 2 cot a; 


2(1 + cot a;) 2 


3 

(8.5556) 


2.5727 


14.25% 


2.9356 


2.14% 


Scarf II 
(hyperbolic) 


a = 1 
A = 1 
B = 1 


tanh a; 
+secha: 


— sech 2 a; 
+3secha;tanhai 


1.0811 


1.19257 


10.31% 


1.02981 


4.74 % 


|-rPT1PT?ll 1 r 7(^f\ 

VJCllCL dllZj t Li 

Poschl- 
Teller 


a = 1 
A = 1 
B = 3 


cothr 
— 3cschr 


2 V1U 

+ llcsch 2 r 
— 9cschr cothr 


1.33772 


1.2790 


4.39% 


1.2576 


5.99% 


Scarf I 
(trigonometric) 


a = 1 
A = 3 
B = 1 


3 tan a; 
— sec a; 


-i-vs 

+7 sec 2 a; 
—5 sec a; tan a; 


3.0505 
(10.0505) 


2.689 


11.8% 


3.1285 


2.56% 



Table 2: Comparison of the exact ground state energies of shape invariant 
potentials with results from the WKB and Friedrich- Trost method. The percent 
errors are also shown. The SWKB results are not shown, since they are always 
exact for the potentials under consideration. 



